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I. INTRODUCTION 



Graphene, a two-dimensional single layer of graphite, was first fabricated in 2004 by 
Novoselov et. al}. This has provided an unique opportunity for experimental observation 
of electronic properties of graphene which has attracted theoretical attention for several 
decades-^. The importance of graphene lies not only in providing the first realization of 
Dirac physics in condensed matter systems but also in providing a way of realization of 
several devices in nanometer scale. In this article, we are going to concern ourselves mainly 
on the first of these two aspects of graphene. 

In graphene, the energy bands touch the Fermi energy at six discrete points at the edges 
of the hexagonal Brillouin zone. Out of these six Fermi points, only two are inequivalent; 
they are commonly referred to as K and K' points^. The quasiparticle excitations about 
these K and K' points obey linear Dirac-like energy dispersion. The presence of such 
Dirac-like quasiparticles is expected to lead to a number of unusual electronic properties in 
graphene including relativistic quantum Hall effect with unusual structure of Hall plateaus^. 
Recently, experimental observation of the unusual plateau structure of the Hall conductivity 
has confirmed this theoretical prediction^. Further, as suggested in Ref. iQ, the presence of 
such quasiparticles in graphene provides us with an experimental test bed for Klein paradox.- 
These and several other properties of graphene has been covered extensively in several review 
articles^""—. In the current article, we are going to focus on the effect of the Dirac nature 
of graphene quasiparticles on two separate aspects. The first of these involves transport of 
superconducting graphene junctions while the second involves Kondo effect and scanning 
tunneling spectra of graphene. 

It is well known that the existence Dirac-like quasiparticles affects tunneling conduc- 
tance of a normal metal-superconductor (NS) interface of graphene^i. Graphene is not a 
natural superconductor. However, superconductivity can be induced in a graphene layer in 
the presence of a superconducting electrode near it via proximity effeclr^"— or by possible 
intercalation with dopant moleculesW. It has been recently predictedrA that a graphene 
NS junction, due to the Dirac-like energy spectrum of its quasiparticles, can exhibit spec- 
ular Andreev reflection in contrast to the usual retro reflection observed in conventional 
NS junctionsi^ii^. Such specular Andreev reflection process leads to qualitatively different 
tunneling conductance curves compared to conventional NS junctionsii. The effect of the 
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presence of a thin barrier region of thickness c? — )■ created by applying a large gate voltage 
Vq ^ oo (such that VqcI is 



also been studied in Ref. 
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inite) between the normal and the superconducting region has 
It has been shown that in this thin barrier limit, in contrast to 
all normal metal-barrier-superconductor (NBS) junctions studied so far, the tunneling con- 
ductance of a graphene NBS junction is an oscillatory function of the dimensionless barrier 
strength x = yod/{fiVF), where vp denotes the Fermi velocity of graphene, with periodicity 
TT. Further, it has also been demonstrated that the tunneling conductance reaches its max- 
ima of 2Go for X = (^ + l/2)7r, where n is an integer. The latter result was also interpreted in 
terms of transmission resonance property of the Dirac-Bogoliubov quasiparticles^. However, 
no such studies have been undertaken for NBS junctions with barriers of arbitrary thickness 
d and barrier potential Vq. As we shall discuss in details in Sec. Ill A[ the analysis of Ref. 



171 and calculate the tunneling conductance of a graphene NBS junction with a barrier of 
thickness d and with an arbitrary voltage Vq applied across the barrier region can also be 
extended to thick barrier junctions^^. It can be shown that the oscillatory behavior of the 
tunneling conductance is not a property of the thin barrier limit, but persists for arbitrary 
barrier width d and applied gate voltage Vq, as long as <^ ^, where ^ is the coherence 
length of the superconductor. Further, the periodicity and amplitude of these oscillations 
deviate from their values in the thin barrier limit and becomes a function of the applied 
voltage Vq. 

The study of Josephson effect in graphene for tunnel SBS junctions also presents some 
unconventional features due to the presence of the Dirac quasiparticles. In this review, we 
shall concentrate on SBS junctions with barrier thickness d <^ ^ where ^ is the superconduct- 
ing coherence length, and width L which has an applied gate voltage Vq across the barrier 
region^. The central property of such junctions on which we shall mainly focus on is that 
in complete contrast to the conventional Josephson tunnel junctions studied so far— i^, the 
Josephson current in graphene SBS tunnel junctions is an oscillatory function of both the 
barrier thickness d and the applied gate voltage Vq. In the thin barrier limit, where the bar- 
rier region can be characterized by an effective dimensionless barrier strength x = V^d/hvp 
[vp being the Fermi velocity of electrons in graphene), the Josephson current becomes an 
oscillatory function of x with period n^^. In this limit, the oscillatory behavior of Josephson 
current can be understood as a consequence of transmission resonance phenomenon of Dirac- 
Bogoliubov-de Gennes (DBdG) quasiparticles in graphene. The Josephson current reaches 



the Kulik-Omelyanchuk limit^^ for x = ('^ being an integer), but, unlike conventional 
junctions, never reaches the Ambegaokar Baratoff limit^ for large x- This analysis is done 
in Sec. ITTBl 

Another extremely interesting phenomenon in conventional metal systems is the Kondo 
effect which occurs in the presence of dilute concentration of localized quantum spins coupled 
to the spin-degenerate Fermi sea of metal electrons^i. The impurity spin-electron interaction 
then results in perfect or partial screening of the impurity spin as one approaches zero 
temperature. It also results in a sharp 'Kondo Resonance' in electron spectral functions. 
Recent developments in quantum dots and nano devices have given new ways in which 
various theoretical results in Kondo physics, which are not easily testable otherwise, can 
be tested and confirmed experiment ally2S. Most of the early studies in Kondo effect were 
carried on for conventional metallic systems with constant density of states (DOS) at the 
Fermi surface^^. Some studies on Kondo effect in possible flux phasea^I, nodal quasiparticles 
in d-wave superconductors^^, Luttinger liquida^S^ and hexagonal Kondo lattice^, for which 
the DOS of the associated Fermions vanishes as some power law at the Fermi surface, has 
also been undertaken. Recently, there has been a interest in study of the physics of magnetic 
impurities in graphene.— One of the purpose of this article is to articulate a part of this 
recent progress in Sec. IIV A[ 

Scanning tunneling microscopes (STM) are extremely useful probes for studying prop- 
erties of two or quasi-two dimensional materials^i^. Studying electronic properties of a 
sample with STM typically involves measurement of the tunneling conductance G{y) for a 
given applied voltage V . The tunneling conductances measured in these experiments have 
also been studied theoretically for conventional metallic systems and are known to exhibit 
Fano resonances at zero bias voltage in the presence of impurities^*^. The application 
and utility of this experimental technique, with superconducting STM tips, has also been 
discussed in the literature for conventional systems^''. However, tunneling spectroscopy of 
graphene using superconducting STM tips remains to be studied both experimentally and 
theoretically. In Sec. lIVBt we shall elaborate the progress on the STM response of doped 
graphene and discuss some of it's unconventional features. For undoped graphene with 
Fermi energy Ep = 0, the derivative of the STM tunneling conductance (G) with respect 
to the applied voltage (dC/dV) reflects the density of states (DOS) of the STM tip (pt), 
i.e., dG/dV ~ +{—)Pt for V > (<)0. By tuning Ep, one can interpolate between this 
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unconventional pt ~ ±dG/dV and the conventional pt ^ G (seen for Ep ^ eV) behaviors. 
Further, for superconducting STM tips with energy gap Aq, G{dG/dV) displays a cusp 
(discontinuity) at eV = —Ep — Aq as a signature of the Dirac point which should be exper- 
imentally observable in graphene with small Ep where the regime eV > Ep can be easily 
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have seen that 



accessed. For impurity doped graphene with large Ep, experiments in Ref. 
the tunneling conductance, as measured by a metallic STM tip, depends qualitatively on 
the position of the impurity in the graphene matrix. For impurity atoms atop the hexagon 
center, the zero-bias tunneling conductance shows a peak; for those atop a graphene site, it 
shows a dip. We provide a detailed discussion of this phenomenon and point out that its 
origin lies in conservation/breaking of pseudospin symmetry of the Dirac quasiparticles by 
the impurity. 



The organization of the rest of the review is as follows. We give a generic description of 
the graphene NBS and SBS tunnel junctions which is described by the Dirac-Bogoliubov-de 
Gennes (DBdG) equations. In section III Al we review the theory of tunneling conductance 
of a graphene NBS junction with a barrier of thickness d and with an arbitrary voltage Vq 
applied across the barrier region. The results obtained are then compared and contrasted 
with that of a thin barrier and zero barrier junction. In section III Bl we study Josephson 
current for a general SBS junction barrier region of thickness d and potential Vq. We also 
discuss the thin barrier limit to understand the oscillatory behavior in terms of transmission 
resonance of DBdG particles. Finally we study some possible experimental realizations of 
the above mentioned junctions to probe the oscillations in section ilTTl In section. IIV Al we 
discuss the unconventional Kondo effect in graphene. We describe the large N analysis for 
a generic spin S local moment coupled to Dirac electrons in graphene. The analysis gives 
rise to a finite critical Kondo coupling strength which can be tuned by the application of an 
external gate voltage and is particular to graphene. We also discuss the possible realization 
of the non- Fermi liquid ground states via the multichannel Kondo effect. In section IIVBI 
we discuss the STM response of graphene. We discuss the tunneling current through the 
STM tip within linear-response theory using a superconding tip to probe an undoped sample 
and a metallic tip with constant density of states (DOS) to probe an impurity present in 
the sample. We conclude with a general discussion on the unconventional tunneling, STM 
properties and the behaviour of magnetic impurities in graphene in section |Vl 
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II. TRANSPORT PROPERTIES OF SUPERCONDUCTING JUNCTIONS 



An understanding of the transport properties across different superconducting junctions 
of graphene throws substantial hght on the electronic properties. A generic description 
of the junctions to study the transport properties is as follows. A local potential barrier 
of width d is implemented on the graphene sheet occupying the xy plane by either using 
the electric field effect or local chemical doping^i^!^ A s-wave pairing A(r) is induced in 
graphene via proximity effeclr^i^. For NBS (SBS) region I is normal (superconducting) 
region occupying x < d for all y as shown schematically in Fig. [H The region II modeled 
by a barrier potential Vq, extends from x = d to x = Q while the superconducting region 
occupies X > (marked as region III in Fig ([T]). For calculations we shall assume that the 
barrier region has sharp edges on both sides. This condition requires that d <^ X = 2?: /kp , 
where kp and A are Fermi wave- vector and wavelength for graphene, and can be realistically 
created in experimentsS. Also, the interface is smooth and impurity free on the scale of the 
superconducting coherence length ^ = hvp/'^o, where Aq is the amplitude of the induced 
superconducting order parameter. For both the junctions (NBS and SBS) the induced 
pair potential can be modeled (with appropriate boundary condition for the two different 
junctions) as: 

A(r) = Aoexp(#) (1) 

(j) is the phase. These junctions can then be described by the Dirac-Bogoliubov-de Gennes 
(DBdG) equations: 

:^-Ep + Uiv) A(r) \ , J,, ... 

Ipa = Elpa (2) 

A*(r) EF-U(v)-naj 

Here, tpa = (ipAayi^Bayi^Aay ~4'Ba) 4 compoucut wavcfunctious for the electron and 

hole spinors, the index a denote K or K' for electron/holes near K and K' points, a takes 
values K'{K) for a = K{K'), Ep denote the Fermi energy. A and B denote the two 
inequivalent sites in the hexagonal lattice of graphene, and the Hamiltonian Ha is given by 

Ha = -ihvp {(Txdx + sgn{a)aydy) . (3) 

In Eq. 131 vp denotes the Fermi velocity of the quasiparticles in graphene and sgn(a) takes 
values ± for a = K[K'). The potential U{y) gives the relative shift of Fermi energies in the 
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barrier and superconducting regions and is modeled as: 



U{r) = Vo9{-x)e{x + d) 



(4) 



Eq. |2] can be solved in a straightforward manner to yield the wavef unction ip in the normal, 
insulating and the superconducting regions, taking into account both Andreev and normal 
reflection processes. These wavefunctions satisfy the appropriate boundary conditions at the 
interfaces of the junctions. Note however that these boundary conditions, in contrast their 
counterparts in standard junction interfaces, do not impose any constraint on derivative 
of the wavefunctions at the boundary. The tunneling conductance and Josephson current 
across the junctions can then be calculated using appropriate expressions. These are found 
to have novel oscillatory behavior in complete contrast to their standard counterparts as 
will be described in the subsequent sections. 

A. NBS junction 

The pair-potential for the NBS junction is modeled as: 



6{x) is the Heaviside step function. Eq. |2] can be solved in a straightforward manner to 
yield the wavefunction in the normal, insulating and the superconducting regions. In the 
normal region, for electron and holes traveling the ±x direction with a transverse momentum 
ky = q and energy e, the (unnormalized) wavefunctions are given by 




(5) 



(l,±e±^0,0) exp[ti±knx + qy)], 




(6) 



where the wave- vector kn{k'^) for the electron (hole) wavefunctions are given by 





and a{a') is the angle of incidence of the electron (hole). 
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Superconducting 
(for SBS junctions) 
A--exp(iOi) 




Superconducting 

{for NBS/SBS junctions) 
A~exp(i02) 



FIG. 1. A schematic sketch of the graphene superconducting junction. Region I denotes the 
normal (superconducting) region for NBS(SBS) junction. A potential barrier Vq of width d is 
created with the application of the external gate voltage. The region III is the superconducting 
region for both the NBS and SBS junctions. Superconductivity, as discussed in the text is induced 
in the grpahene layer by proximity effect. 

In the barrier region, one can similarly obtain 

= (1, ±e^^', 0, 0) exp {±hx + qy)] , 
V'^^ = (O, 0, 1, Te±^^') exp [i {±klx + qy)] , (8) 

for electron and holes moving along ±x. Here the angle of incidence of the electron (hole) 
9{9') and the wavevector kb{k'^) are given by is 

sin [^(^0] = fivpq/ [e + (-){Ep - Vb)] , 



Note that Eq. [8] ceases to be the solution of the Dirac equation (Eq. [2]) when Ep = Vq and 
e = 0. For these parameter values, Eq. [2] in the barrier region becomes 'Ha'^B = which 
do not have purely oscillatory solutions. For the rest of the calculation, we shall restrict 
ourselves to the regime Vq > Ep. 

In the superconducting region, the DBdG quasiparticles are mixtures of electron and 
holes. Consequently, the wavefunctions of the DBdG quasiparticles moving along ±x with 
transverse momenta q and energy e, for [Uq + Ep) ^ Aq, e, has the form 

= [wf, wf, w^, w^) eW±'=^-+''2^)+'^-l (10) 

where 

= ±exp(±i7), = exp[-i(0i T P)], 

^ = ±exp[±z(T0i + /3 + 7)], (11) 

where 7 is the angle of incidence for the quasiparticles. Here the wavevector kg and the 
localization length can be expressed as a function of the energy e and the transverse 
momenta q as 
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ks = \J[{Uo + Ef) /hvpY -q 

[(t/o + i?i.)Aosin(/3)]' ^ ^ 



where /3 is given by 



/3 = cos"^ (e/Ao) if |e| < Aq, 
= -i cosh-^ (e/Ao) if |e| > Aq. (13) 

Note that for |e| > Aq, k becomes imaginary and the quasiparticles can propagate in the 
bulk of the superconductor. 

Next we note that for the Andreev process to take place, the angles ^, 9' and a' must all 
be less than 90°. This sets the limit of maximum angle of incidence a. Using Eqns. Eland 
[9l one finds that the critical angle of incidence is 

ac = a^^^e{Vo - 2Ep) + af^e{2EF - Vq) 
= arcsin [|e - Ep\ / {e + Ep)] , 

= arcsin [|e - \Ep - Vo\\ / {e + Ep)] . (14) 



FIG. 2. A schematic sketch of normal reflection (r), Andreev reflection (r^) and transmission 
processes (t and t') at a graphene NBS junction. Note that in this schematic picture, we have 
chosen to denote a retro Andreev reflection for iflustration purpose. In practice, as discussed 
in the text, takes into account possibihties of both retro and specular Andreev reflections. The 
electron and hole wavefunctions inside the barrier region is not sketched to avoid clutter. 



Note that in the thin or zero barrier limits treated in Refs. 171 and lUJ, ac = a^c^ for all 



parameter regimes. 

Let us now consider a electron-like quasiparticle incident on the barrier from the normal 
side with an energy e and transverse momentum q. The basic process of ordinary and 
Andreev reflection that can take place at the interface is schematically sketched in Fig. |2J 
As noted in Ref. Ill iii contrast to conventional NBS junction, graphene junctions allow for 
both retro and specular Andreev reflections. The former dominates when e, Aq <^ Ep so 
that a = —a' (Eq. [6]) while that latter prevails when Ep <^ e, Aq with a = a'. Note that in 
Fig. m we have chosen to denote a retro Andreev reflection for illustration purposes. In 
practice, va includes both retro and specular Andreev reflections. In what follows, we shall 
denote the total probability amplitude of Andreev reflection as which takes into account 
possibilities of both retro and specular Andreev reflections. 

The wave functions in the normal, insulating and superconducting regions, taking into 
account both Andreev and normal reflection processes, can then be written as^ 

= P^'it + #b" + ^^b'^ + (15) 

where r and are the amplitudes of normal and Andreev reflections respectively, t and 
t' are the amplitudes of electron-like and hole-like quasiparticles in the superconducting 
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region and p, q, m and n are the amplitudes of electron and holes in the barrier. These 
wavefunctions must satisfy the appropriate boundary conditions: 

^7vU=-d = ^BU=-d, ^b|x=0 = ^5|x=0. (16) 

These boundary conditions yield eight linear homogeneous equations for the coefficients r, 
r^, t, t' , p, q, m, and n. 

After some straightforward but cumbersome algebra, we find that 

r = (17) 

AT = [e^" cosihd + e)-i sm{khd)] 

-p[cos{kbd -6) -i e*" sm.{k},d)], (18) 
V = [e"^" cos{hd + e)+i sm{hd)] 

+p [cos{hd -9)+ ie-*" sm{hd)] , (19) 

= cos(g)[re-/^ + e-/^] (['°'^^^^ - ^) - 

+re'''"'^[cos{hd -6)+ ie"*" sin(A;ferf)]) , (20) 

t = n', (21) 

t(r + l)e*'="'='cos(^')e"''^ 
cos{k'f^d - 6*') - ie-*"' sin(A;^(i) 

where the parameters T and p can be expressed in terms of 7, /3, 9, 9', a, and a' (Eqs. El [HI 

[ini and [ED as 

_re«(7-^) _|_ g-«(7-/3) 

r = (24) 

_ e"*"' cos(A;;rf + 9') - i sin(fc^ci) 
^ ~ cos(fc^ci - 9') - ze-*"' sin(fc;rf) ' ^ ^ 

The tunneling conductance of the NBS junction can now be expressed in terms of r and ta 
G{eV) r f. , ,2 , , ,2 cos(a') 



GoieV) Jo 



r(l-\rf + \rAr-^)cosia)da, 
Jo V cos(a) J 

(26) 



where Gq = Ae^N{eV)/h is the ballistic conductance of metallic graphene, eV denotes the 
bias voltage, and A^(e) = {Ep + e)w/ {nhvp) denotes the number of available channels for 
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a graphene sample of width w. For eV <^ Ep^ Gq is a constant. Eq. [26] can be evaluated 
numerically to yield the tunneling conductance of the NBS junction for arbitrary parameter 
values. We note at the outset, that G = when ac = 0. This occurs in two situations. 
First, when eV = Ep and Vq > 2Ep so that = ai'^ vanishes. For this situation to arise, 
Ep + f/o > A > Ep which means that Uq has to be finite. Second, = etc = when 
eV = and Ep = Vq, so that the zero-bias conductance vanishes when the barrier potential 
matches the Fermi energy of the normal side — 

We now make contact with the results of the thin barrier limit. We note that since 
there are no condition on the derivatives of wavefunctions in graphene NBS junctions, the 
standard delta function potential approximation for thin barrier^® can not be taken the 
outset, but has to be taken at the end of the calculation. This limit is defined as rf/A — )■ 
and Vo/Ep — )■ oo such that the dimensionless barrier strength 

X = Vod/hvp = 2n (^) (^) (27) 



^Ep J 

remains finite. In this limit, as can be seen from Eqs. El [9] and [TOl 6, 6', knd, k'^d — > and 
k^d, k'^d — )• X so that the expressions for F, p and rj (Eq. [25]) 

ptb ^ e"^'^ - T]^^ ^ e"^"' cos(x) -isin(x) 

gi7 _|_ ^tb ' cos(x) — ie~*"' sin(x) ' 

where the superscript "tb" denotes thin barrier. Using the above-mentioned relations, we 
also obtain 



r 



„tb 



cos(x) (e*" — p*^) — zsin(x) (l — p*'^e*") 
cos(x) (e-*" + p"^) + i sin(x) (1 + p^^e"*") ' 
'tb cos(x) (l + r^^) — i sin(x) (e*° — r^^e""^) 



t 



^tb ^ r^'tb 



^ cos(x) - «e-*" sm(x) 



Eqs. [28] and [29] are precisely the result obtained in Ref. [l7| for the tunneling conductance of 
a thin graphene NBS junction. The result obtained for a zero barrier in Ref. can be now 
easily obtained from Eqs. [28] and [29] by substituting x = in these equations. 
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1. Qualitative Discussions 



In this section, we shall analyze the formulae for tunneling conductance obtained in above 
section. First we aim to obtain a qualitative understanding of the behavior of the tunneling 
conductance for finite barrier strength. To this end, we note from Eq. [22] that the maxima 
of the tunneling conductance must occur where |rp is minimum. In fact, if |rp = for all 
transverse momenta, the tunneling conductance reaches its value 2Gq. Therefore we shall 
first try to analyze the expression of r (Eq. [T7|) for subgap voltages and when the Fermi 
surfaces of the normal and superconducting sides are aligned with each other (f/g = 0). In 
this case, we need Aq <^ Ep- So for subgap tunneling conductance, we have e < Aq ^ Ep. 
In this limit, a ~ —a' ~ 7 (Eqs. and [TUj) . fc^ ~ /c^, and 6 ~ —9' (Eq. [H]). Using these, one 
can write 

6*°" cos{kbd — 6) — i sm{khd) 



V 
P 



cos{kf,d + 6) — ie*" sin{kf,d) 
T] cos(a — (3) + i sin(/3) 



(30) 
(31) 



cos(a + l3) + irj sin(/3) 

Substituting Eq. [3T] in the expression of A/", we find that the numerator of the reflection 
amplitude r becomes (Eqs. [T7] and [TS]) 



e' 
Do 



— 4sin(a) sin(/3) cos{kbd — 



X — i cos(a) sm{kbd) 

(cos(M -0) + cos(M + ^))/2 
2 [cos{kbd + 6) - cos{kbd - 6)] 
cos(a — (3) {cos(a) + [cos{kbd 



X 



(32) 



+ cos(M + 0)] /2} + sm{kBd) sin(/3) 

Dq = cos{kbd + 9) cos{a + (3) + sin{kbd) sin(/3) 

+ie^°' [cos{kbd — 0) sin(/3) — sm{kbd) cos(a + /?)] . 

(33) 

From the expression of M (Eq. [521) . we note the following features. First, for normal 
incidence (a = 0) where 6 = 6' = 0, and hence r (Eq. [T71) vanishes. Thus the barrier is 
refiectionless for quasiparticles which incident normally on the barrier for arbitrary barrier 
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FIG. 3. Plot of zero-bias tunneling conductance for Uq = and Aq = O.OlEp as a function of gate 
voltage Vo and barrier thickness d. Note that the oscillatory behavior of the tunneling conductance 
persists for the entire range of Vq and d. 

thickness d and strength of the apphed voltage Vq. This is a manifestation of Klein paradox 
for Dirac-Bogoliubov quasiparticles^. However, this feature is not manifested in tunneling 
conductance G ( Eq. which receives contribution from all angles of incidence. Second, 
apart from the above-mentioned cases, r never vanishes for all angles of incidence a and 
arbitrary eV < Aq unless 9 = 9'. Thus the subgap tunneling conductance is not expected 
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eV/A, 



FIG. 4. Plot of tunneling conductance of a NBS junction graphene as a function of bias voltage for 
different effective barrier strengths for Uq = and Aq = O.OlEp. Note that the curves for x = 
(black line) and x = 7r(pink circles) coincide reflecting vr periodicity. 

to reach a maximum value of 2Go as long as the thin barrier limit is not satisfied. However, 
in practice, for barriers with Vq > AEp, the difference between 6 and 6' turns out to be 
small for all q < kp {< 0.25 for q < kp and eV = 0) so that the contribution to (Eq. |32|) 
from the terms ~ {cos{ki,d + 6) — cos{ki)d — 6)) becomes negligible. Thus |rp can become 
quite small for special values of Vq for all g < A;^ so that the maximum value of tunneling 
conductance can reach close to 2Go- Third, for large Vq, for which the contribution of terms 
~ {cos{kiid + 6) — cos{khd — 9)) becomes negligible, and hence r becomes very small when 
the applied voltage matches the gap edge i.e. sin(/3) = (Eq. [32]). Thus the tunneling 
conductance curves approaches close to its maximum value 2Gq and becomes independent 
of the gate voltage Vq at the gap edge eV = Aq for Aq ^ Ep, as is also seen for conventional 
NBS junctions^^. Fourth, in the thin barrier limit, {Vq/Ep — )■ oo and d/X 0), 6 ^ and 
kbd — ^ X) so that the contribution of the terms ~ {cos{kbd + 9) — cos{khd — 9)) in Eq. |32] 
vanishes and one gets 



.rtb ^ 2 sin(a) [sin(x + P) - sin(x - f3)] 

ntb 

X [- cos(x) + i sin(x) cos(a)] , 



tb 



cos(x) cos(a + f3) + sin(x) sin(/3) + ie* 
X [cos(x) sin(/3) — sin(x) cos(q: + (3)] . 



(34) 
(35) 
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FIG. 5. Plot of periodicity Xperiod of oscillations of tunneling conductance as a function of applied 
gate voltage Vq for Uq = and Aq = 0.01£'_f. Note that the periodicity approaches vr as the voltage 
increases since the junction approaches the thin barrier limit. 

A/"*^ and hence r**^ (Eq. [29]) vanishes at x = [n + l/2)7r which yields the transmission 
resonance condition for NBS junctions in graphene and is given in Fig@^. 

Fifth, as can seen from Eqs. ITTl and l22l both |rp and Ir^P are periodic functions of Vq 
and d since both and 6 depend on Vq. Thus the oscillatory behavior of subgap tunneling 
conductance as a function of applied gate voltage Vq or barrier thickness d is a general 
feature of graphene NBS junctions with d <^ C,- However, unlike the thin barrier limit, for 
an arbitrary NBS junction, k^d = x\/ {Ef/Vq — 1)^ + h^v'^pQ^ /Vq ^ x-, and Q ^ ^. Thus 
the period of oscillations of |rp and will depend on Vq and should deviate from their 
universal value vr in the thin barrier limits^. Finally, we note from Eqs. [TTJ |26] and [3l] that 
in the thin barrier limit (and therefore for large Vq), the amplitude of oscillations of the zero- 
bias conductance for a fixed Vq, defined as [G'max(eV^ = 0; Vq) —G^^J^eV = 0; Vo)]/G'o! which 
depends on the difference of \r{x = {n + l/2)7r)p and \r{x = n7r)p becomes independent of 
X or the applied gate voltage Vq. 

2. Numerical Results 

The above-mentioned discussion is corroborated by numerical evaluation of the tunneling 
conductance as shown in Figs. El [5], [6] and [71 From Fig. [3l which plots zero-bias tunneling 
conductance G{eV = 0) as a function of Vq and d, we find that G{eV = 0) is an oscillatory 
function of both Vq and d and reaches close to its maximum value of 2Gq throughout the 
plotted range of Vq and d. Further, as seen from Fig. |5l the periodicity of these oscillations 
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FIG. 6. Plot of the amplitude [Gmax(eF = 0;Vo) - G^ir^ieV = 0;Vo)]/Go = (G„,ax - G^i^)/Go 
of zero-bias tunneling conductance as a function of the applied gate voltage Vq for Uq = and 
Aq = O.OlEp. Note that G reaches 2Go for Vq > 4Ef where the amplitude become independent 
of the applied gate voltage as in the thin barrier limit and vanishes for Vo/Ep = 1 as discussed in 
the text. 

becomes a function of Vq. To measure the periodicity of these oscillations, the tunneling 
conductance is plotted for a fixed Vq as a function of d. The periodicity of the conductance 
(^period is uoted dowu from these plots and Xperiod = Vodperiod/^vp is computed. Fig. |5]clearly 
shows that Xperiod deviate significantly from their thin barrier value vr for low enough Vq 
and diverges at Vq — )■ Ep^ . Fig. |6] shows the amplitude of oscillations of zero-bias conduc- 
tance as a function of Vq. We note that maximum of the zero-bias tunneling conductance 
Graaxi^V = 0) reachcs close to 2Go for Vq > Vqc — 4:Ep. For V > Vqc, the amplitude 
becomes independent of the applied voltage as in the thin barrier limit, as shown in Fig. |6l 
For Vq — )■ Ep, = ttc — )• 0, so that G{eV = 0) — i- and hence the amplitude vanishes. 
Finally, in Fig. [TJ we plot the tunneling conductance G as a function of the applied bias- 
voltage eV and applied gate voltage Vq for d = 0.4A. We find that, as expected from Eq. [Ml 
G reaches close to 2Gq at the gap edge for all Vq > 6Ep. Also, as in the thin barrier limit, 
the oscillation amplitudes for the subgap tunneling conductance is maximum at zero-bias 
and shrinks to zero at the gap edge eV = Aq, where the tunneling conductance become 
independent of the gate voltage. 

Next, we consider the case Uq 7^ 0, so that Aq — Ep <^ {Ep + Uq). In this regime, 
there is a large mismatch of Fermi surfaces on the normal and superconducting sides. Such 
a mismatch is well-known to act as an effective barrier for NBS junctions. Consequently, 
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FIG. 7. Plot of tunneling conductance as a function of the bias- voltage eV and gate voltage Vq for 
d = 0.4A and Aq = O.OlEp. Note that for large Vq, the tunneling conductance at eV = Aq is close 
to 2Go and becomes independent of Vq (see text for discussion). 



additional barrier created by the gate voltage becomes irrelevant, and we expect the tunnel- 
ing conductance to become independent of the applied gate voltage Vq. Also note that at 
eV = Ep, ttc = (Eq. [T^ . Hence there is no Andreev reflection and consequently Gq van- 
ishes for all values of the applied gate voltage for this bias voltage. Our results in this hmit, 
coincides with those of Ref. Finally in Fig. [9l we show the dependence of amplitude of 
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FIG. 8. Plot of tunneling conductance as a function of the bias- voltage eV and the gate voltage 
Vo for d = 0.4A, Aq = 2Ep and Uq = 2bEp. As discussed in the text, the tunneling conductance 
is virtually independent of the applied gate voltage Vq due to the presence of a large C/q. Note 
that maximum angle of incidence for which Andreev reflection can take place vanishes at eV = Ep 
leading to vanishing of G at this bias voltage. 
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FIG. 9. Plot of amplitude of oscillation (Gmax — G'min)/G'o of zero-bias tunneling conductance as 
a function of Uq/Ef for Vq = QEp and Aq = Q.QlEp- The oscillation amplitudes always decay 
monotonically with increasing Uq independent of Vq. 

oscillation of zero-bias tunneling conductance on Uq for the applied bias voltages Vq = QEp 
and Aq = O.OlEp. As expected, the oscillation amplitude with decreases monotonically with 
increasing Uq. We have verified that this feature is independent of the applied gate voltage 
Vq as long as Vq > Vqc- 

19 



B. SBS junction 



For this junction, the region I in Fig [T] is a superconducting region and the pair-potential 
can be given as 

A(r) = Ao [exp{i(t)2)9ix) + exp(i0i)^(x + d)] (36) 

where Aq is the amphtude and 0i(2) are the phases of the induced superconducting order 
parameters in regions I (II) as shown in Fig{T], and 6 is the Heaviside step function. Solving 
Eq. [21 the wavefunctions in the superconducting and the barriers regions are obtained. In 
region I, the wavefunctions for the DBdG quasiparticles moving along ±x direction with a 
transverse momentum ky = q = 27m/ L (for integer n) and energy e, are given by^ 

^± = {uf,ut, uf, u^) e^i(^''^^+iy)+--] (37) 

where 

^ = ± exp(±i7), = expH(0i T /?)], 

u-^ Ui 

^ = ±exp[±z(T0i + /3 + 7)], (38) 

and X]i=i 4 I — 2k is the normalization condition for the wavefunction for d <C k"^, where 
= (/j,fi?)^/cs/ [£'i?Ao sin(/3)] is the localization length. Here kg = \J [Ep/hvpf' — (f, 7, 
the angle of incidence for the quasiparticles, is given by sin(7) = hvpq/Ep, and (3 is given 
by 

13 = cos"^ (e/Ao) if |e| < Aq, 
= -I cosh"^ (e/Ao) if |e| > Aq, (39) 

Note that for |e| > Aq, k becomes imaginary and the quasiparticles can propagate in the 
bulk of the superconductor. The wavefunctions in region II (x > ) can also be obtained 
in a similar manner 

ijjfj = (y± y±, vf, v^) eW±^»^+92')-«^], (40) 
where J2i=i 4 l^jP = 2k and the coefficients Vi are given by 

4 = ±exp(±i7), 4 = expH(02±/3)], 



4 = ±exp[±2(T02-/3 + 7)], (41) 



20 



The wavefunctions for electrons and holes moving along ±x in the barrier region is given 

by 

= (1, ±e±'^ 0, 0) exp [i {±hx + qy)] /V2d, 
= (^0, 0, 1, ±e^^^') exp [t {±k'f^x + qy)] /V2d. 

Here the angle of incidence of the electron(hole) 0{6') and is given by: 

hvpq 



sin [9{9')] 



e+i-)iEF-Vo) 



hiK)-^|r-^i^¥^] -f (42) 



F - Vq) ^ ^ 

To compute the Josephson current in the SBS junction, the energy dispersion of the 
subgap Andreev bound states are found which are localized with localization length at 
the barrier— The energy dispersion e„ (corresponding to the subgap state characterized by 
the quantum number n) of these states depends on the phase difference = 02 — 0i between 
the superconductors. The Josephson current / across the junction at a temperature Tq is 
given hy^^ 

mx,n) = ^j: f: ^/(e„), (43) 

n q=~kp 

where f{x) = l/(e^/(^-s7b) _|_ jg ^^ie Fermi distribution function and is the Boltzman 
constant^ 

To obtain these subgap Andreev bound states, boundary conditions at the barrier are 
imposed. The wavefunctions in the superconducting and barrier regions are constructed 
using Eqs. [371 140] and |42] as 

where 01(02) and 61(62) are the amplitudes of right and left moving DBdG quasiparticles in 
region I(II) and p{q) and r(s) are the amplitudes of right(left) moving electron and holes 
respectively in the barrier. These wavefunctions satisfy the boundary conditions: 

^/U=-d = ^ijU=-d, ^bU=o = ^//U=o- (45) 
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Substituting Eqs. [STJ SQl HH andHHin Eq. HSl we get eight linear homogeneous equations 
for the coefficients ai=i,2, &i=i,2, P, r, and s, so that the condition for non-zero solutions 
of these coefficients can be obtained as 



A'sm{2/3) + B' cos(2/3) + C = 



(46) 



where A! , B', andC are given by 



A' = cos{k'fjd) cos(7) cos{9') sm{kbd) (sin(7) sin(6') — 1) 
+ cos{kbd) cos(7) cos(6') sin(A;^(i) 
+ - cos{kbd) cos{9) sin(27) sin(6'') sin(A;[,(i) 

B' = sin(fc[,(i) sm{kbd) [ — 1 + sin(6') sin(7) 

— sin(^') sin(7) + sin(6') sin(6'') sin^(7)] 

— cos{khd) cos{k'f^d) cos^(7) cos(6') cos(6'') 

C' = cos^(7) cos(6') cos(6'') cos(0) — sm(kbd) sm{k'bd) 
X [sin(^) sin(6'') - sin^(7) 



Note that in general the coefficients A', B', and C depends on e through kb, k^, 6 and 9' which 
makes it impossible to find an analytical solution for Eq. HHl However, for subgap states in 
graphene SBS junctions, e < Aq ^ Ep. Further, for short tunnel barrier we have | Vq — -E^l ^ 
Ep- In this regime, as can be seen from Eqs. |42l A', B', and C become independent of e 



since /cf, ~ A;[, ~ fci = a/[(-E'f — Vo)/hvF]'^ — and 6 ~ —6' 6i = sin ^ [hvpg/iEp — Vq)] 
so that the e dependence of kb, fc^, and 6' can be neglected. In this regime one finds that 
A', B', C -^A,B,C where 



+ sin(7) (sin(^) - sin(e'))] 



(47) 



A 







B 



C 



— sin^lkid) [1 — sin(7) sin(6'i)]^ 

— cos'^{kid) cos^(7) cos^(6'i) 
sm^{kid) [sin(7) — sin(6'i)]^ 
+ cos^(7) cos^(^i) cos(0) 



(48) 
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FIG. 10. Plot of Josephson current / as a function of phase difference (j) and the appHed gate 
voltage Vq for ksTo = O.OIAq and d = 0.5A showing oscillatory behavior of I/Iq as a function of 
the applied gate voltage. 



The dispersion of the Andreev subgap states can now be obtained from Eqs. |16] and [391 
There are two Andreev subgap states with energies e± = ±e where 



e = Ao^/l/2-C/2B (49) 

Using Eq. |33l one can now obtain the expression for the Josephson current 

I {(f), Vo, d, To) = /o5f(0, Vo, d, Tq), 

cos^(7) cos^(6'i) sin(0) 



Vo, d, To) = d'j 



Bt/^o 
X tanh(e/2A;ijro) 



(50) 



where /q = eAoEFL/2h'^nvF and we have replaced J2q ~^ EpL/ (2TxhvF) J^^J^^dj cos{'y) as 
appropriate for wide junctions^^. 

The dispersion of the Andreev subgap states and the Josephson current in graphene SBS 
junctions, in complete contrast to their conventional counterparts^Si^il^^ is found to be an 
oscillatory function of the applied gate voltage Vq and the barrier thickness d. This statement 
can be most easily checked by plotting the Josephson current J as a function of the phase 
difference and the applied gate voltage Vq for a representative barrier thickness d = 0.5A 
and temperature ksTo = O.OIAq, as done in Fig. [101 In Fig. [HI we plot the critical current 
of these junctions /c(Vo, d, Tq) = Max[/(0, Vq, d, Tq)] as a function of the applied gate voltage 
Vq and barrier thickness d for low temperature /csTq = O.OIAq. The critical current of these 
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FIG. 11. Plot of Ic/Io vs the applied gate voltage Vb and the junction thickness d for To = O.OIAq. 

graphene SBS junctions is an oscillatory function of both Vq and d. This behavior is to be 
contrasted with those of conventional junctions where the critical current is a monotonically 
decreasing function of both applied bias voltage Vq and junction thickness cPSiSiiiS. 

We analyze a few other properties of these oscillations. To find the amplitude of oscil- 
lation, we compute Ic as a function of Vq (for a representative value of d = 0.3A), note the 
maximum (I™*^^) and minimum (/™™) values of Ic, and calculate the amplitude J™^^ — /™™. 
The procedure is repeated for several temperatures Tq and the result is plotted in Fig. [T2] 
which shows that the amplitude of oscillations decreases monotonically as a function of 
temperature. Next, we discuss the period of oscillation of the critical current. To obtain 
the period, the critical current Ic as a function of barrier width d for the fixed applied gate 
voltage Vo is computed and (iperiod is noted down . Then Xperiod = Vorfperiod/^^^F is computed 
and Xperiod as a function of Vq for /c^^o = O.OIAq is plotted as shown in Fig. [131 It is found 
that Xperiod decreases with Vq and approaches an universal value vr for large Vq > 20Ep. 
This property, as we shall see in the next section, can be understood by analysis of graphene 
SBS junctions in the thin barrier limit (Vq — )■ oo and d ^ such that x = Vod/hvp re- 
mains finitei^) and is a direct consequence of transmission resonance phenomenon of DBdG 
quasiparticles in superconducting graphene. 

1. Thin barrier limit 

In the limit of thin barrier, where Vq ^ oo and — )■ such that x = V^d/hvp remains 
finite, 6*1 —t- and kid — )■ x- From Eqs. Il8]and|l9l we find that in this hmit, the dispersion 
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FIG. 12. Plot of the temperature dependence of the amphtude of oscihations of Ic (given by 
[I^^^{d) — /™°((i)]//o) for d = 0.3A. The amphtude is measured by noting the maximum and 
minimum values of the critical current by varying Vq for a fixed d. 

of the Andreev bound states becomes 



et{q, <\>- x) = ±Ao VI - ^(7, x) sin2(0/2), 
cos^(7) 



T(7,X) 



1 — cos2(x) sin^(7) 



(51) 
(52) 



where the superscript 'tb' denote thin barrier hmit. The Josephson current / can be obtained 
substituting Eq. |52] in Eq. |43l In the limit of wide junctions, one gets 



/*'(0,X,To) 



9'\<P,X,To) 



W"(0,x,To 



7r/2 



tt/2 



d'-f 



T(7, x) cos{'y) Aschemat sin(0) 
v/l-T(7,x) sin^(0/2) 



X tanh {e+/2kBTo) 



(53) 



We find that the Josephson current in graphene SBS junctions is a vr periodic oscillatory 
function of the effective barrier strength x the thin barrier limit. Further we observe 
that the transmission probability of the DBdG quasiparticles in a thin SBS junction is given 
by r(7, x) which is also the transmission probability of a Dirac quasiparticle through a 
square potential barrier as noted in Ref. [Q. Note that the transmission becomes unity 
for normal incidence (7 = 0) and when x = nvr. The former condition is a manifestation 
of the Klein paradox for DBdG quasiparticles^. However, this property is not reflected 
in the Josephson current which receives contribution from quasiparticles approaching the 
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FIG. 13. Plot of 

Xpcriod of the critical current Ic as a function of Vq. Note that Xperiod approaches 
TT as we approach the thin barrier limit. 

junction at all angles of incidence. The latter condition {x = nn) represents transmis- 
sion resonance condition of the DBdG quasiparticles. Thus the barrier becomes completely 
transparent to the approaching quasiparticles when x = nvr and in this limit the Josephson 
current reduces to its value for conventional tunnel junctions in the Kulik-Omelyanchuk 
limit: 1*^(0, nvr, Tq) = 4/o sin(0/2)Sgn(cos(0/2)) tanh (Aq |cos(0/2)| /2A;bTo)22. This yields 
the critical Josephson current lf^{x = nn) = AIq for fc^To ^ Aq. Note, however, that in 
contrast to conventional junctions T(7, x) can not be made arbitrarily small for all 7 by 
increasing x- Hence If^ never reaches the Ambegaokar-Baratoff limit of conventional tunnel 
junctions^^. Instead, if^ix) becomes a vr periodic oscillatory function of x- The amplitude 
of these oscillations decreases monotonically with temperature. 

Finally, the product I^^Rn which is routinely used to characterize Josephson junctions^Si^i 
is computed, where -R^v is the normal state resistance of the junction. For graphene SBS 
junctions Rn corresponds to the resistance of a Dirac quasiparticle as it moves across a 
normal metal-barrier-normal metal junction. For short and wide junctions discussed here, 
it is given by Rn = Rq/si{x) where Rq = n'^vpf? I {e^EpL) and is given by^*^ 

Slix) = c?7^(7,x)-cos(7). (54) 

J~n/2 

Note that si{x) and hence Rn is an oscillatory function of x with minimum 0.5-Ro at x = 
and maximum 0.75i?o at x = -|- l/2)7r. The product I^^Rn, for thin SBS junctions is 
given by 

II'^Rn = {7:^,/2e)gt.^{x, T)/s,{x), (55) 
where Sfmaxlx) is the maximum value of g*'^{(p, x)- Note that I^^Rn is independent of Ep and 
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FIG. 14. Plot of II^Rn as a function of x- Ic -^N is an oscillatory bounded function of x 
never reaches its value (7rAo/2e) for conventional junctions in the Ambegaokar-Baratoff limit. 

hence survives in the limit Ep — 0^^. For /cbTq <^ Aq, fi'max(^'^) = 4 and si{nn) = 2, so that 
Il^RN\x=mT = 7rAo/e which coincides with Kulik-Omelyanchuk limit for conventional tunnel 
junctions2^i^. However, in contrast to the conventional junction, I^^Rn for graphene SBS 
junctions do not monotonically decrease to the Ambegaokar-Baratoff limit^^i^ of 7rAo/2e ~ 
1.57Ao/e as x is increased, but demonstrates vr periodic oscillatory behavior and remains 
bounded between the values vrAo/e at x = "^tt and 2.27Ao/e at x = + l/2)7r, as shown 
in Fig. [H 



III. EXPERIMENTS 

Superconductivity has recently been experimentally realized in graphene^. In the pro- 
posed experiment to observe the oscillatory behavior in the tunneling conductance and 
Josephson current, one needs to realize these junctions in graphene. The local barrier can 
be fabricated using methods of Refs. [5I 41I. The easiest experimentally achievable regime 
corresponds to Aq ^ Ep with aligned Fermi surfaces for the normal and superconducting 
regions. We suggest measurement of tunneling conductance curves at zero-bias {eV = 0) 
in this regime. Our prediction is that the zero-bias conductance will show an oscillatory 
behavior with the bias voltage. In graphene, typical Fermi energy can be Ep < 40meV and 
the Fermi-wavelength is A > 100nn>^>^'^'^. Effective barrier strengths of < 80meV- and 
barrier widths of ~ 10 — 50 nm therefore specifies the range of experimentally feasible 
junctions^i^i^. Consequently for experimental junctions, the ratio Vq/Ep can be arbitrarily 
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large within these parameter ranges by fixing Vq and lowering Ep. Experimentally, one can 
set 5 < Ep < 20meV so that the conditions Aq ^ Ep Vo/Ep 1 is easily satisfied for 
realistic Aq ~ O.SmeV and Vq = 200meV. This sets the approximate range Vo/Ep > 10 for 
the experiments. Note that since the period (amplitude) of oscillations increases (decreases) 
as Vo/Ep — )■ 1, it is preferable to have sufficiently large values of V^/Ep for experimental 
detection of these oscillations. 

To check the oscillatory behavior of the zero-bias tunneling conductance, it would be 
necessary to change Vq in small steps 6Vq. For barriers of a fixed width, for example with 
values of d/X = 0.3, it will be enough to change Vq in steps of approximately 20 — 30meV, 
which should be experimentally feasible. 

We note that for the above-mentioned range of Vo/Ep, the experimental junctions shall 
not always be in the thin barrier limit. For example, as is clear from Fig. [5l the periodicity 
of oscillations Xperiod of the zero-bias tunneling conductance of such junctions shall be a 
function of Vq and shall differ from vr. This justifies the theoretical study of NBS junctions 
in graphene which are away from the thin barrier limit. 

Apart from the above-mentioned experiments, it should also be possible to measure the 
tunneling conductance as a function of the applied bias voltage eV/ Aq for different applied 
gate voltages Vq. Such measurements can be directly compared with Fig. [61 Finally, it 
might be also possible to create a relative bias Uo between the Fermi surfaces in the normal 
and superconducting side and compare the dependence of oscillation amplitudes of zero-bias 
tunneling conductance on Uo with the theoretical result shown in Fig. [HI 

We also suggest measuring DC Josephson current in these junctions as a function of the 
applied voltage Vq. Such experiments for conventional Josephson junctions are well-known^. 
Further SNS junctions in graphene has also been recently been experimentally created^i^. 
To observe the oscillatory behavior of the Josephson current, alike the procedure to measure 
the tunneling conductance, it would be necessary to change Vb in small steps 6Vo. For 
barriers with fixed d/X = 0.3 and Vo/Ep = 10, this would require changing Vq in steps of 
approximately 30meV which is experimentally feasible. The Joule heating in such junctions, 
proportional to I^Rn, should also show measurable oscillatory behavior as a function of Vq. 
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IV. KONDO EFFECT AND STM SPECTRA 



A. Kondo effect in Graphene 

In this section, we shall present a large N analysis for a generic local moment coupled 
to Dirac electrons in graphene to show that Kondo effect in graphene is unconventional can 
be tuned by gate voltage. We demonstrate the presence of a finite critical Kondo coupling 
strength in neutral graphene. We point out that local moments in graphene can lead to non 
Fermi-liquid ground state via multi channel Kondo effect. 

The crucial requirement for occurrence of Kondo effect is that the embedded impurities 
should retain their magnetic moment in the presence of conduction of electrons of graphene. 
We will not quantitatively address the problem of local moment formation in the presence 
of Dirac sea of electrons in graphene in the present paper. We expect that large band 
width and small linearly vanishing density of states at the fermi level in graphene should 
make survival of impurity magnetic moment easier than in the conventional 3D metallic 
matrix. A qualitative estimate of the resultant Kondo coupling can be easily made consid- 
ering hybridization of electrons in vr band in graphene with d orbitals of transition metals. 
Typical hopping matrix elements for electrons in vr band is t ~ 2eV and effective Hubbard 
U in transition metals is 8eV. So the Kondo exchange J ~ At'^/U, estimated via standard 
Schrieffer-Wolf transformation, can be as large as 2 eV which is close to one of the largest 
J ^ 2.5 eV for Mn in Zn. Therefore it is customary to use Kondo Hamiltonian^^ to study 
the effect in Graphene. 

1. Large N analysis 

Our analysis begins with the Hamiltonian for non-interacting Dirac electron in graphene. 
In the presence of a gate voltage V, the Hamiltonian can be expressed in terms of electron 
annihilation operators ^^(b)^ at sublattice A{B) and Dirac point s = K,K with spin a =t, i 

as 





eV hvpikx — isgrL{s)ky) 
hvF{k^ + isgn{s)ky) eV 




(56) 
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where sgn(s) = 1(— 1) for s = K{K'), vp is the Fermi velocity of graphene, and all repeated 
indices are summed over. In Eq. |56]and rest of the analysis, we shall use an upper momentum 
cutoff kc = A/{hvF), where A ~ 2eV corresponds to energy up to which the linear Dirac 
dispersion is valid, for all momenta integrals. 

Eq. inSlcan be easily diagonalized to obtain the eigenvalues and eigenfunctions of the Dirac 
electrons: E± = eV ± hvpk where k = {kx,ky) = {k,6) denote momenta in graphene and 



l/-\/2 (1, ±exp (isgn(s)6')). Following Ref. 



271 . we now introduce the ^ fields, 



which represents low energy excitations with energies E±, and write 

nM = Y.'^AiU = i/v^(ajk) +c.(k)), 

nM = exp{^e)/V2{eU^) - CM)- (57) 

In what follows, we shall consider a single impurity to be centered around x = 0. Thus 
to obtain an expression for the coupling term between the local moment and the conduction 
electrons, we shall need to obtain an expression for \I'(x = 0) = ^^(O). To this end, we 
expand the ^ fields in angular momentum channels ^+a(k) = X]m=-oo where we 

have written k = {k,6). After some straightforward algebra, one obtains 

= 7f /'^ ^ {cT^'^'ik) - c-T^'^'ik)) , 

Note that \1/_b(0) receives contribution from m = ±1 channel while for \1/a(0), the m = 
channel contributes. The Kondo coupling of the electrons with the impurity spin is given 
by 

= ^ E E E E <ioK,^uo)s% (59) 

c s=l 1=1 a,(3=l a=l 

where g is the effective Kondo coupling for energy scales up to the cutoff A, S denotes 
the spin at the impurity site, r are the generators of the SU(Ai"c) spin group, and we have 
now generalized the fermions, in the spirit of large N analysis, to have A^^ flavors (valley 
indices) Nf colors (sublattice indices) and Nc spin. For realistic systems Nf = Nc = Ns = 2. 
Here we have chosen Kondo coupling g to be independent of sublattice and valley indices. 
This is not a necessary assumption. However, we shall avoid extension of our analysis to 
flavor and/or color dependent coupling term for simplicity. Also, the Dirac nature of the 
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graphene conduction electrons necessitates the Kondo Hamiltonian to mix m = ±1 and 
m = channels (Eqs. |58] and |59|) . This is in complete contrast to the conventional Kondo 
systems where the Kondo coupling involves only m = angular momentum channel. 

The kinetic energy of the Dirac electrons can also be expressed in terms of the ^ fields: 

^ E (60) 



m=—co s,a 



Typically such a term involves all angular momenta channels. For our purpose here, it will 
be enough to consider the contribution from electrons in the m = 0, ±1 channels which 
contribute to scattering from the impurity (Eqs. [58] and [59]) . To make further analytical 
progress, we now unfold the range of momenta k from (0, oo) to (— oo, oo) by defining the 

fields C^,(2)a 

cW(^) = +(-)v^^-^''^^^^^^(l^l)' ^<0' (61) 

so that one can express the \l/ fields as '^A{B)a(^) ~ I^oo ^ V^'^i(2)a(k) ■ In terms of the 
c^2)a fields, the kinetic energy (in the m = 0, ±1 channels) and the Kondo terms in the 
Hamiltonian can therefore be written as 



Ho = r dk/{2Tr)Ek4lc^,^ 

J ~hc 

Hk = g/iSn^kl) f r ^/\k\VW\dkdk' 



x(c?J(fc)r,",c?,(A:')5"), (62) 

where E]. = eV + hvpk and summation over all repeated indices are assumed. 

We follow standard procedure^ of representing the local spin by SU(A''c)Fermionic fields 
fa and write the partition function of the system in terms of the / and c fields 

Z = j VcVc^VfVpVe e-^/^ S = Sq + Si + S2 



I3h 



I3h 



So= I dr I dk/{2n){c',^{k,T)Go'clik 



Si = J/{A7T^N,kl) dr / ^/\k\^/\F\dkdk' 

Jo J -kc J -kc 

X r) clp{k', r)fl{T)r:;sfs{T)_ 

S2= dr [{flir) [hd, + e(r)] /,(r)) - e(r)g] , (63) 
Jo 
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FIG. 15. Sketch of the critical Kondo coupling Jc{q,T) as a function of temperature for several 
applied voltages q = eV/A. The Kondo phase exists for J > Jc- 

where Gq^ = Kdr + Ek is the propagator for c fields, J = gNc/2 is the renormalized Kondo 
couphng, we have imposed the impurity site occupancy constraint J2a fafa = Q using a 
Lagrange multipher field €{t). 
We now use the identity^ 

= NcSasSp^ - (64) 

and decouple 5*1 using a Hubbard-Stratonovitch field 0f. In the large N^. limit one has 

5 = 5*0 + 52 + 5*3 + S'4, where 




S, = Kkl/J dT<Prir)<Pt{T). (65) 
Jo 

Note that at the saddle point level (0^*) ~ ^Xla ^llf^^ that a non-zero value of 0f indicates 
the Kondo phase. In what follows, we are going to look for the static saddle point solution 
with 4>^i{t) = 00 and e(r) = etr^. In this case, it is easy to integrate out the c and / fields, 
and obtain an effective action in terms of 0o and eo and one gets S' = S5 + Sq with 

^5 = -/3^iV,Tr [In {tfiLOn - eo - N^Njcj^lG'^itoo, V^)0o)] , 

^6 = I3h [NsN.Nfkl |0o|' I J - eoQ) , (66) 
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where Tr denotes Matsubara frequency sum as well as trace over all matrices and the Fermion 
Green function Gq(?p„, q) = Gq is given by— 

^ 27r(n!L)2 ^^^" ~ - ' (67) 

where, in the last line we have switched to dimensionless variables p„ = hwn/A and q = eV/A. 

One can now obtain the saddle point equations from Eq. [66] which are given by 6S'/6(f)Q = 
and 6S'/6eo = 0. Using Eqs. [66]and|671 one gets (after continuing to real frequencies and 
for T = 0) 

1/J = -A/{7rhvFkiy / dpGoip -V- AoGo/2)-\ 

g/iVe = l/(27r) y dpv{p~v- /\oGo/2)-\ (68) 

where we have defined the dimensionless variable Aq = N fNs\(t)o\'^ / {irfi^vy), p = hu/A, 
Go = 2n{hvF)'^G'Q/ A, v = to/A > 0, and have used the energy cutoff A for all frequency 
integrals. At the critical value of the coupling strength, putting z/ = and Aq = 0, we 
finally obtain the expression for Jdq, T) 

J,{q, T) = J,(0) [1 - 2q In (l/g^) In (^A)] ' (69) 

where the temperature ksT is the infrared cutoff, Jc(0) = {nhypk^Y / A = vr^A is the critical 
coupling in the absence of the gate voltage, and we have omitted all subleading non-divergent 
term which are not important for our purpose. For \^ = = g, we thus have, analogous to 
the Kondo effect in flux phase systems^l, a finite critical Kondo coupling Jc(0) = vr^A ~ 20eV 
which is a consequence of vanishing density of states at the Fermi energy for Dirac electrons 
in graphene. Of course, the mean-field theory overestimates Jc- A quantitatively accurate 
estimate of Jc requires a more sophisticated analysis which we have not attempted here. 



2. Results and Discussions 



The presence of a gate voltage leads to a Fermi surface and consequently Jc{q, T) — )• as 
T — 7- 0. For a given experimental coupling J < Jc(0) and temperature T, one can tune the 
gate voltage to enter a Kondo phase. Fig. [151 which shows a plot of Jc{q, T) as a function of 
T for several gate voltages q illustrates this point. The temperature T*(g) below which the 
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system enters the Kondo phase for a physical couphng J can be obtained using Jc(g, T*) = J 
which yields 



ksT* = Aexp [(1 - J,(0)/J)/(2gln[l/g2])] (70) 

For a typical J ~ 2eV and voltage eV ~ 0.5eV, T* ~ 35K— We stress that even with 
overestimated Jc, physically reasonable J leads to experimentally achievable T* for a wide 
range of experimentally tunable gate voltages. 

We now discuss the possible ground state in the Kondo phase qualitatively. In the absence 



of the gate vo^ 
noted in Ref. 



tage a finite Jc implies that the ground state will be non-Fermi liquid as also 
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for flux phase systems. In view of the large Jc estimated above, it might 
be hard to realize such a state in undoped graphene. However, in the presence of the gate 
voltage, if the impurity atom generates a spin half moment and the Kondo coupling is 
independent of the valley(fiavor) index, we shall have a realization of two-channel Kondo 
effect in graphene owing to the valley degeneracy of the Dirac electrons. This would again 
lead to overscreening and thus a non Fermi-liquid like ground state^^. The study of details 
of such a ground state necessitates an analysis beyond our large mean-field theory. To 
our knowledge, such an analysis has not been undertaken for Kondo systems with angular 
momentum mixing. In this work, we shall be content with pointing out the possibility of 
such a multichannel Kondo effect in graphene and leave a more detailed analysis as an open 
problem for future work. 

Next, we discuss experimental observability of the Kondo phenomena in graphene. The 
main problem in this respect is creation of local moment in graphene. There are several 
routes to solving this problem, i) Substitution of a carbon atom by a transition metal atom. 
This might in principle frustrate the strong sp^ bonding and thus locally disturb the integrity 
of graphene atomic net. However, nature has found imaginative ways of incorporating 
transition metal atoms in p-vr bonded planar molecular systems such as porphyrin^. Similar 
transition metal atom incorporation in extended graphene, with the help of suitable bridging 
atoms, might be possible, ii) One can try chemisorption of transition metal atoms such as Fe 
on graphene surface through sp-d hybridization in a similar way as in intercalated graphite^, 
iii) It might be possible to chemically bond molecules or free radicals with magnetic moment 
on graphene surface as recently done with cobalt pthalocyanene (CoPc) molecule on AU(lll) 
surface^. This might result in a strong coupling between graphene and impurity atom 
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leading to high Kondo temperatures as seen for CoPc on AU(lll) surface {Tk — 280K). 
iv) Recently ferromagnetic cobalt atom clusters with sub nano-meter size, deposited on 
carbon nanotube, have exhibited Kondo resonance^. Similar clusters deposition in graphene 
might be a good candidate for realization of Kondo systems in graphene. v) From quantum 
chemistry arguments, a carbon vacancy, or substitution of a carbon atom by a boron or 
nitrogen might lead to a spin-half local moment formation. In particular, it has been shown 
that generation of local defects by proton irradiation can create local moments in graphite^. 
Similar irradiation technique may also work for graphene. 

For spin one local moments and in the presence of sufficiently large voltage and low 
temperature, one can have a conventional Kondo effect in graphene. The Kondo temperature 
for this can be easily estimated using ksTx ~ Dexp( — 1/pJ) where the band cutoff D ~ 
lOeV, J ~ 2 — 3eV and DOS per site in graphene p ^ 1/20 per eV. This yield Tx — 6 — 150K. 
The estimated value of Tk has rather large variation due to exponential dependence on J. 
However, we note that Kondo effect due to Cobalt nano-particle in graphitic systems such 
as carbon nanotube leads to a high Tx ~ 50K which means that a large J may not be 
uncommon in these systems. Recently, It has also been shown that the Kondo effect can 
be controlled by orbital degrees of freedom™. A symmetry class of orbitals in a magnetic 
adatoms with inner shell in graphene leads to a distinct quantum critical points, where 
Tk oc ( J — JcY^^ near the critical coupling J^. 

Finally, we note that recent experiments have shown a striking conductance changes in 
carbon nanotubes and graphene, to the extent of being able to detect single paramagnetic 
spin-half NO2 molecule^. This has been ascribed to conductance increase arising from hole 
doping (one electron transfer from graphene to NO2). Although Kondo effect can also lead 
to conductance changes, in view of the fact that a similar effect has been also seen for 
diamagnetic NH3 molecules, the physics in these experiments is likely to be that of charge 
transfer and not local moment formation. 

B. STM spectra of Graphene 

The experimental situation for STM measurement is schematically represented in Fig. 
[T6l The STM tip is placed atop the impurity and the tunneling current X is measured as 
a function of applied bias voltage V. The possible positions of the impurity is shown in 
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FIG. 16. Schematic experimental setup with the right panel showing two possible positions (atop 
hexagon center and atop a B site) of the impurity. The numbers denote nearest neighbor A and 
B sublattice sites to the impurity. ai(2) = +(— )\/3/2£ + 3/2y [lattice spacing set to unity] are 
graphene lattice vectors. The choice of coordinate center (0,0) are shown for each case 

the right panel of Fig. Such a situation can be modeled by the well-known Anderson 
Hamiltonian^. Here we incorporate the low-energy Dirac quasiparticles of graphene in this 
Hamiltonian which is given by 



H = HG + H, + Ht + HGd + Hat + (71) 

Hg= I i^l\k) [hvF{T,a,K + (Tyky) - EpI] V^f (fc) 
Jk 

Hd= ^ eddlds + Un^n^ 



Ht = Y.[Y. + (^ot^t^L^ + h-c)] (72) 

E / fc(^)<.(^)4 + h.c.) (73) 



a=A,B • 



Had 

Hdt= Yl {WXsdl + h.c). (74) 
HGt= [ (uUkKsik)tl + h.c.) (75) 

Here Hq is the Dirac Hamiltonian for the graphene which are described by the two 
component annihilation operator ips{k) = (c^^(A;), c'^gik)) belonging to the valley P = K,K' 
and spin s =t, i, I is the identity matrix, r and a denotes Pauli matrices in valley and 
pseudospin spaces, vp is the Fermi velocity, and = X]/3=_ft: k' ^s=t4 / (Cp"' denotes 
the impurity atom Hamiltonian with an on-site energy and U is the strength of on- 

36 



site Hubbard interaction. Ht is the Hamiltonian for the superconducting (Aq 7^ 0) or 
metaUic (Aq = 0) tip electrons with on-site energy en, where v signifies all quantum numbers 
(except spin) associated with the tip electrons. The operators dg and ti,s are the annihilation 
operators for the impurity and the tip electrons. The Hamiltonians HqcLi Hct-, and Hcn 
describe interaction between the graphene and the impurity electrons, the graphene and 
the STM tip electrons, and the impurity and the STM tip electrons respectively. The 
corresponding interactions parameters V^{k)^ U^.^iji), and are taken to be independent 
of valley and spin indices of graphene electrons but may depend on their sublattice index or 
pseudospin. In defining the Hamiltonian Eq. [TH we omit inter-valley scattering of electrons 
by the impurity. This is usually justified if the impurity radius is larger than \K — K'\^^ so 
that the inter-valley scattering is suppressed. In the present case, for impurity atom atop 
the hexagon center, the impurity scattering potential respects the pseudospin symmetry 
(since it does not distinguish between A 01 B sites) and hence can not fiip pseudospin of 
graphene electrons. Thus a graphene electron in momentum state ki around K valley can 
only be scattered to around momentum states ^2 = —ki in K' valley. This constraint further 
reduces the phase space for inter-valley scattering. This phase space constraint makes the 
inter-valley scattering small for all impurities atop hexagon center. For impurities atop 
graphene sites, our analysis is applicable for large impurity size where inter- valley scattering 
is negligible. 



1. Tunneling Current 



The tunneling current for the present model is given by 



X{t) = e{dNt/dt) = ie{[H,Nt])/h 



(76) 



where N = J2i^s A^s^ys is the number operator for the tip electrons. These commutators 
receive contribution from Hdt and Hat in Eqs- EH and [75] and can be evaluated by a straight- 
forward generalization of method outlined in Ref. Isol to the case of superconducting tips. A 
standard calculation yields 



X = Xn 



x|[/°P + 



duj[f{uj - eV) - f{uj)]pt{uj - eV) pg{uj 
\B{co)\' \q{u)\' - 1 + 2Re[q{uMuj)- 



ImSrf(cc;) 



(77) 
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where Iq = 2e(l + C^)/h, pai^) and pt{e) are the graphene and STM tip electron DOS 
respectively, C, = \U^\/\Uy^\ = is the ratio of coupling of the impurity to the 

electrons in B and A sites of graphene with f/^ = f/° and = V^, x(e) = [e — — 
ReS(i(e)]/lmSfj(e),/(e) = 1/(1 + exp[e/T]) (/c^ = 1) is the Fermi function, and ^d{^) is the 
impurity advanced self-energy in the absence of the tip. Here B{e) = V^U^l2{t) and g(e) is 
given by 

^^'^ v^) ^^^^ 

where, 

h{e)={l + e)Y.^r[Re{g{e,k)}] (79) 

k 

h{e)={l + e)Y.^r[lm{g{e,m (80) 

k 

Here, we have neglected the energy dependence of the coupling functions assuming small 
applied voltages. Tr denotes trace over Pauli matrices in pseudospin, valley and spin spaces, 
and Q is the Green function for the graphene electrons: 

(e + Ef)I — hVpiTzCTx 

hj- -\- Cyhyj /Ol \ 



A simple calculation yield; 



/,(e) = -4(1 + e')(e + E^) In |1 - AV(e + Epf] /A^ 

/2(e) = 4(1 + e)rr\e + Ep\e{A - e - Ep)/k\ (82) 

where A is the ultraviolet momentum cutoff and 9 is the Heaviside step function. Usually, in 
graphene, A is taken to be the energy at which the graphene bands start bending rendering 
the low-energy Dirac theory inapplicable and can be estimated to be 1 — 2eV— . 



2. Results and Discussions 

Here we are going to analyse the tunneling conductance G = as measured by STM. 
In the absence of impurities, the contribution to the conductance comes from the first 
term of Eq. [771 For s-wave superconducting tips, one finds that the tunneling conductance 
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FIG. 17. (Color online) Plot of the tunneling conductance G and its derivative dG/dV as a function 
of the applied bias voltage eF/Ao = —p for r = 0,2,6 (red solid, blue dashed and black dotted 
lines) respectively. See text for details. 



{G{V) = dX/dV) for > and at T = is given by (with r = Ep/^o, p = -eV/Ao 



G = Go Aft{p)\r\+ / Sgn{z - p + r)Aft{z)dz 



(83) 



dC _ eCo 
dV~^ 



M(p) -Kip)\r \ - 26'(p - r)Mt{p - r) 



(84) 



where Go 

\x' ' 



87rV|f/op(l + e)potpo/h, Pg = Po\r - p\, pt{r) = potM(r), Aftix) 



AN- 



I — 1), and Sgn(x) denote the signum function. For graphene with Ep = 
r = 0, dG/dV ~ Sgn{V)Aft{—V), i.e., the tip DOS is given by the derivative of the tun- 
nehng conductance. For large Ep away from the Dirac point, the first term of G becomes 
large and reflects the tip DOS. In between these extremes, when Ep ~ eV, neither G nor 
dG/dV reflects the DOS. In this region, the signature of the Dirac point appears through a 
cusp (discontinuity) in G (dG/dV) at eV = —Ep — Aq arising from the contribution of the 
second (third) term in Eq. |83] (Eq. These features, shown in Fig. [T71 distinguishes such 
graphene STM spectra with their conventional counterparts^. 

Next, we turn to the case of impurity doped graphene and consider a metallic tip with 
constant DOS. The contribution to the tunneling conductance from the impurity (after 
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FIG. 18. (Color online) Plot of Gimp as a function of V for |14/^'^/C/'^| = 0.05 (right; impurity atop 
a site) and 2 (left; impurity atop hexagon center) for Ep/A = 0.3,0.1, andO (black solid, blue 
dashed and red dotted lines respectively). Plot parameters are 5U = V° = 0.05A, W° = 0.0005A, 
and ed = 0. 



subtracting the graphene background) at T = (Eq. |82]) is 



^ _ ^' \BiV)\' \qiV)\' - 1 + 2Re[qiV)]xiV) 

'^^-'^-'^°ImS,(V) A[1 + X^(V^)] ' ^ ^ 

where Gg = 2e^pQt^/h. Such tunnehng conductances are known to have peak/antiresonace/dip 
feature at zero bias for |g| ^ 1/ ~ 1/ ^ 1^^. In conventional metals, Eq. 1771 can be used to 
compute the STM current by taking as a fixed parameter independent of the position of 
the impurity. However, the situation in graphene necessitates a closer attention to which 
is proportional to the probability amplitude of the Dirac quasiparticles in graphene to hop 
to the tip. The strength of f/° can be estimated using the well-known Bardeen tunneling 
formula;^: 

t/° ~ y rfV {ct>l{z)d,^G{r, z) - ^5,(f, z)d,ct>,{z)) ~ ^G(ro, ^o) (86) 

where the last similarity is obtained by a careful evaluation of the surface integral J (Pr 
over a surface between the graphene and the tip parallel to the graphene sheet^, (ro, zq) 
is the coordinate of the tip center-^^, (l)y{z) is tip electron wavefunction, and the wavefunc- 
tion graphene electrons \[^G(r, z) around K{K') valley, can be written, within tight-binding 
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approximation, as^ 

^G(r, = ^ E e^'^^^^'^-'''^-^-'^ [fir- Rf) 

+e+^-)is^^{r-Rf)\f{z). (87) 

— * 

Here 9^ = aTctan{ky/kx), Sk is the Fermi wave- vector as measured from the Dirac points with 
\6k\ <^ \K{K') \ for all Ep, (p{r) are localized pz orbital wavef unctions, is a normalization 
constant, f{z) is a decaying function of z with decay length set by work function of graphene, 
and Rf^^^ = nai +77102(0.2 ~y) with integers n, m denote coordinates of the graphene lattice 
sites (Fig. [T6l) ^°. When the impurity and the STM tip is atop the center of the hexagon, 
pseudospin symmetry necessitates f{fo — Rf'^) to be identical for all neighboring A and B 
sublattice points 1..6 surrounding the impurity (Fig. [12]). Consequently, the sum over lattice 
vectors Rf in Eq. [87] reduces to a sum over the phase factors exp(i[{i^'(i^'') + 6k} ■ Rf]) for 
these lattice points. It is easy to check that this sum vanishes for both Dirac points (when 
\6k\ = 0). Thus the only contribution to "^Gi'^o^^o) comes from the second and further 
neighbor sites for which the amplitude of localized wavefunctions ip{fo — R^ ) are small. 
For finite Ep, {6k 7^ 0) there is a finite but small contribution (0(|5A;|/|K|)) to "^ci^OyZo) 
from the nearest neighbor sites. Thus \l/G(ro, zq) and hence is drastically reduced when 
the impurity is atop the hexagon center. In this case, we expect f/° ^ and hence 
|g| ^ 1 (Eq. 178]) leading to a peaked spectra for all Ep. In contrast, for the impurity 
atom atop a site, there is no such symmetry induced cancellation and ipci'f^OjZo) receives 
maximal contribution from the nearest graphene site directly below the tip. Thus we expect 
|f/°| ^ (since it is easier for the tip electrons to tunnel to delocalized graphene band 
than to a localized impurity level) leading to g ^ /1//2 — — In |1 — A^/ {eV + Ep)'^\/tt. For 
large |e\^ + i?ir| and impurity atop a site, g < 1 leading to a dip or an antiresonance in Gimp 
which is qualitatively distinct from the peaked spectra for impurity atop the hexagon center. 
As Ep — )■ 0, g diverges logarithmically for small eV. However, it can be shown that in this 
regime x shows a stronger linear divergence for eV 7^ which suppresses Gimp- A.t cV — 6,^, 
the divergence of x also becomes logarithmic and we expect a peak of Gimp- Note that 
these effects are independent of and hence of the precise nature of the impurity. Such 
an impurity position dependent peak/dip structure of Gimp has been observed for magnetic 



impurities in Ref. 
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for Ep > eV. 
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To demonstrate this feature, we restrict ourselves to impurities with small Hubbard U 
and compute the self energy of the impurity electrons within a mean-field theory where 
Uria-n^ = U{n^)ns leading to spin-dependent on-site impurity energy = ed + U{n^)^. 
Using Eqs. [7T] and [731 one then obtains the mean-field advanced impurity Green function 
= (w — Co- — T,d{uj))^^ where the impurity self-energy is given by = | + 

i/2) and mean-field self-consistency condition demands = j duo / TT\m.Q^^^ {u) . Following 



Ref. 



571 . we solve these equations to get and ImSd(e) which can be substituted in Eq. [85] 



to obtain Gimp- We note, from Eqs. [85] and [78l that Gimp/Gg depends on the ratios Ep/A, 
V^/A, and W^/U^ which can not be quantitatively determined from the Dirac- Anderson 
model. We therefore treat them as parameters of the theory and compute Gimp for their 
representative values as shown in Fig. [TS] In accordance with earlier discussions, we find 
that for large Ep/A = 0.3, Gimp has qualitatively different features; for the impurity at the 
center of the hexagon, it shows a peak (left panel) while for that atop a site (right panel), it 
shows a dip. The change of Gimp from a dip to a peak via an antiresonance as a function of 
Ep/A when the impurity is atop a site can be seen from right panel of Fig. [18] In contrast, 
the left panel always shows peak spectra. 

Before ending this section, we note that the logarithmic divergence of the Gimp when the 
impurity is reasonably close to the Dirac point is a characteristics of the Dirac physics of the 
low-energy quasiparticles. This feature is therefore also expected to be seen for tunneling 
conductance measured atop an impurity on the surface of a topological insulator. 



V. CONCLUSION 

In this review we have presented a theory for transport properties across superconducting 
junctions of graphene with barriers of thickness do and arbitrary gate voltages Vq applied 
across the barrier region. The oscillatory behaviour of the tunneling conductance as well as 
Josephson current are shown to be robust even for a barrier of finite width. In the thin bar- 
rier limit, such behavior is the manifestation of the transmission resonance of DBdG quasi- 
particles in superconducting graphene. Graphene is an interesting candidate for transport 
applications, in particular for spintronics as it exhibits exhibits remarkably high mobility 
with easily controllable carrier density. Superconducting junctions of graphene has recently 
been realized experimentally^!^. Further experiments in this direction may therefore lead 
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to realization of SBS and NBS junctions discussed in this review and hence may lead to 
verification of some the theoretical results discussed here. 

The effect of localized impurities on the electronic properties of graphene has attracted 
a lot of recent attention. We have studied the effect of presence of localized magnetic 
impurities in graphene which gives rise to the Kondo effect i.e. the dynamic screening of 
the localized moment. The Kondo effect in graphene is unconventional as the effective 
coupling Kondo coupling strength (for weak coupling regime) can be tuned by gate voltage. 
Recent studies^ have found that the Kondo coupling strength can also be controlled by a 
gate voltage in the strong coupling regime. We also discuss scanning tunneling conductance 
spectra phenomenon for both doped and undoped graphene. The position of the impurity 
on or in graphene plays a subtle role and affects the underlying physics of STM spectra 
in doped graphene. For impurity atoms atop the hexagon center, the zero-bias tunneling 
conductance shows a peak; for those atop a graphene site, it shows a dip. This feature is 
a direct consequence of pseudospin symmetry and Dirac nature of graphene quasiparticles. 
A recent scanning tunneling microscopy (STM) experiment^ has demonstrated the ability 
to controllably ionize individual Co adatom on graphene using either a back gate voltage 
or the STM tip bias voltage. This has opened up the possibility of probing the interesting 
electronic phenomena which arises due to the interplay of the impurity with the graphene 
electrons. 

We thank I. Paul and S. Bhattacharjee for collaboration on related topics. K.S thanks 
DST for support through grant SR/S2/CMP-001/2009. 
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